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, In a chiral superconductor with broken time-reversal symmetry a "spontaneous Hall effect" 

' may be observed. We analyze this phenomenon by taking into account the surface properties of 

04 ' a chiral superconductor. We identify two main contributions to the spontaneous Hall effect. One 

, contribution originates from the Bernoulli (or Lorentz) force from spontaneous currents running 

^ ' along the surfaces of the superconductor. The other contribution has a topological origin and is 

I related to the intrinsic angular momentum of Cooper pairs. The latter can be described in terms 

of a Chern-Simons-like term in the low-energy field theory of the superconductor and has some 
similarities with the quantum Hall effect. The spontaneous Hall effect in a chiral superconductor 
is, however, non-universal. Our analysis is based on three approaches to the problem: a self- 
' consistent solution of the Bogoliubov-de Gennes equation, a generalized Ginzburg-Landau theory, 

O ' and a hydrodynamic formulation. All three methods consistently lead to the same conclusion that 

^ , the spontaneous Hall resistance of a two-dimensional superconducting Hall bar is of order h/ {ekpX)^ , 

5-H ' where kp is the Fermi wave vector and A is the London penetration depth; the Hall resistance is 

substantially suppressed from a quantum unit of resistance. Experimental issues in measuring this 
^ . effect are briefly discussed. 
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I. INTRODUCTION 



Unconventional superconductivity appears with a large variety of possible phases displaying many properties that 
are not shared by conventional s-wave superconductors. These phases are characterized by their symmetry properties, 
and in most cases not only the U(l)-gauge symmetry but other symmetries are also spontaneously broken. Interesting 
^ ' physics emerges, in particular, when time-reversal symmetry T is violated in a superconductor. Volovik and Gor'kov 
, have classified such superconducting states into two categories, the so-called "ferromagnetic" and the "antiferromag- 
^ ^ netic" states.EJ They are distinguished by the internal angular momentum of Cooper pairs. In the ferromagnetic state 
, the Cooper pairs possess either a finite orbital or (for non-unitary states) spin moment, while in the antiferromagnetic 
■ state they have no net moments. Examples of these two types of states were recently discussed in connection with 
T-H " high-temperature superconductors; the so-called (1^2 _y2 -\- idxy-'^&ve state represents a ferromagnetic pairing state, 
C3 while the d^^-y^ A- is-wave state is antiferromagnetic. In high-temperature supercouductors a T-violating state is 
most likely to be realized only near surfaces or interfaces at very low temperatures.cTcl Among the heavy fermion 
superconductors there are two well-known systems which have T-violating bulk superconducting phases: UPts and 
Ui_2;Th2,Bei3 (0.017 < x < 0.45). These materials show superconducting double transitions, and T- violation is 
L associated, with the second of the two transitions^ A more recent candidate for T-violating superconductivity is 
Sr2Ru04,u a layered perovskite compound. Experimentally, Sr2Ru04 shows clear features of a strongly correlated 
Q ' quasi-two-dimensional Fermi liquid above the onset temperature of superconductivity Tc =1.5 K. It was suggestedQ 
O ■ that, since this system in the normal state behaves like a two-dimensional analogue to "^He, thepSuaerconducting phase 
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would correspondJxLthe superfluid A-phase, which is a well-known T-violating pairing stateS Ii3 Indeed, the exper- 
imental evidepcdlilEj for unconventional superconductivity in Sr2Ru04 is consistent with the basic order parameter 
^ ■ symmetryJl2Hi3 

d{k) = zA '^^, \ (1.1) 

kp 

a p- wave (spin-triplet) state. Here we have used the standard notation of the d- vector to represent the order parameter 



of the triplet state:la A(A;) — idlk) ■ cra^ , where a* are the Pauli matrices. The rf- vector (IT) pointing along the z 
direction implies that the spin part of the Cooper pair wave function is the spin-triplet state with Sz — 0, i.e., in-plane 
equal-spin pairing (the z direction is along the c-axis of Sr2Ru04). In a system with cylindrical symmetry the orbital 
part of the pair wave function is a state with finite angular momentum along the z-axis, = ±1. Obviously the 



state represented by Eq. (1.1) is "ferromagnetic," and is also called chiral p-wave state. 
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The chirality of this superconducting state is characterized by a topological number N defined b 

oo OG 

1 /" [ „ ^ f 9rri dm 



(1.2) 



The unit vector fn{kx, ky) is 



TTl 

m=- — r, m = {Redz,lmd^,ek) , (1.3) 
m 



where efc — (fc^ + ky)/2m — /i is the kinetic energy measured from the chemical potential fj,. (Throughout this paper 
we will assume the cylindrical symmetry whenever microscopic modeling is necessary, although the inclusion of crystal 
field effects is straightforward. Note that, being topological, the number N can be defined without resorting to the 
cylindrical symmetry.) The topological number N corresponds to the multiplicity of the wrapping of a sphere S'^ in 
the mapping of U {oo} ~ to another S'^ representing the unit vector m. The vect or m does not vanish at any 



fe, allowing us to define the topological number N uniquely. In the chiral p-wave state (1.1) the topological number 
N is either +1 or —1 depending on the chirality of the state or the orbital angular momentum. 

A chiral superconductor has gapiess chiral edge modes at interfaces between states of opposite chiralities (domain 
walls) or between a chiral state and a vacuum (surface), as dictated by the index theorem. Volovik pointed outt£l that 
the topological invariant N determines the number of gapiess edge modes per spin at a boundary between different 
states.E3 This yields 2|A^i — iV2| such modes at the interface of two domains characterized by the topological numbers 
Ni and N2, respectively. At an interface between a chiral p-wave state and fl_.V|acuum we find two modes (one per 
spin). These states can be easily understood in terms of Andreev bound^tatesO^EII Recent point contact experiments 
provide the first indication for these states at the surface of Sr2Ru04.c3 

The presence of gapiess chiral edge modes is reminiscent of quantum Hall fiuids. One can thus expect a chiral 
superconductor to have a Hall effect. This is a spontaneous Hall effect (SHE) in the sense that a transverse voltage 
appears in response to an external current even without an external magnetic field applied to the superconductor The 
Hall conductance is not quantized, however, contrary to a naive analogy. This is because the electric current can be 
carried as a supercurrent by supjepionducting condensate, or, ip-nthjex words, the local charge density is not a conserved 
quantity in a superconductorll^ oil In fact, it has been shownEj'CJ'EJ that the spin and the thermal Hall conductance 
are quantized to the topological number N in appropriate units. Even though it is non- universal, the spontaneous 
(charge) Hall effect is interesting in its own right because it is an effect which can be experimentally measured, in 
principle. Thus we would like to investigate in detail in this paper the origin and magnitude of the effect. To this 
end, it is necessary to understand what one would actually measure in an experiment trying to observe the SHE. 
Computing a diagram for the current-current correlation function in a bulk superconductor does not give a proper 
solution to this problem. On the contrary, we have to understand physics near a surface of a chiral superconductor, 
taking into account issues such as (1) the spatial variation of the order parameter near the surface, (2) the spontaneous 
current running along the surface of the chiral superconductor, and (3) screening effects for both magnetic field and 
electric field in the superconductor. We study the surface properties using three different approaches: a numerical 
self-consistent analysis of the Bogoliubov-de Gennes (BdG) equation, an extended Ginzburg-Landau (GL) theory with 
a scalar potential, and a phenomenological hydrodynamic formulation. We find that there are two basic contributions 
to the SHE. The first is the Bernoulli force due to spontaneous currents which flow near the surface. They are the 
consequence of the chiral surface states, but can also be understood in terms of a transverse order parameter textures. 
Naturally they lead to a Hall response. The second originates from the intrinsic magnetic moment of the Cooper pairs 
in the chiral p-wave state. The former contribution is present in both chiral and non-chiral superconductors, since 
both can have spontaneous surface currents, while the latter contribution exists only in chiral superconductors. We 
therefore conclude that even a non-chiral (antiferromagnetic) T-violating superconductor can have a SHE. 

The organization of this paper is as follows. In Sec. II we introduce a simple microscopic model for the chiral 
p-wave superconductor and solve the BdG equation self-consistently at zero temperature. In Sec. HI we present the 
extended GL theory and derive basic equations. They are used in Sec. IV to analyze the SHE near Tc. In Sec. V 
we present parameters appearing in the extended GL theory obtained from a microscopic weak-coupling analysis. 
The GL theory is the basis to formulate an effective hydrodynamical theory presented in Sec. VI, which allows us to 
interpret and estimate the SHE in a simple way distinguishing the two contributions mentioned above. In Sec. VII 
we also demonstrate the SHE in the non-chiral s + id-wave state by the numerical BdG approach. The results are 
summarized in Sec. VIII. 
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II. SELF-CONSISTENT BOGOLIUBOV-DE GENNES ANALYSIS 



A. Model formulation 



We consider a chiral p-wave superconductor in which energy band and order parameter have no momentum depen- 
dence in the z direction. We may thus treat the system as if it is two-dimensional. The starting point of our analysis 
is a mean- field Hamiltonian for a two-dimensional spin-triplet superconductor with d || 2. The Hamiltonian is given 
by -ffMF = J dx J dy H-mf with the Hamiltonian density 



V 9 



V'|(T-)V'0|(r) + V'|(r)VV'|(T 
1 ^ 

-A{r) -fi-eAoir), 
c J 



T]* (r) • (r) VV'T (r) + i/j^ (r) V-0| (r)] 



(2.1) 
(2.2) 



where ipa- is the annihilation operator of electron with spin a =t, i, .9 is the coupling constant of the attractive 
interaction that is responsible for p-wave pairing (g > 0), and V = {dx,dy). The superconducting order parameter 
should satisfy the (self-consistence) gap equation, obtained from {5 / 5r)*){HMF) = 0, 



17 = {rix.Tlv) 



2fcp 



(^T W[VV'i(r)] - [V^T(r)]V^x(^)) 



We are interested in system s wi th boundaries in which the order parameter r/ depends on r 
system with d given in Eq. (1.1) corresponds to the case r] — A(l, ±i). 

The equations of motion for the field il^a are ihdti'a 
into them, we arrive at the BdG equation: 



(2.3) 

{x,y). The uniform 



[■ipa,HMF]- Substituting tp^ = ue and V'l 



ve 



-iet/h 



ho{r)u„{r) - -^rr—ir] ■ Vu„ + V • (w„J7)] = e„M„(r), 
Zkp 



-ho*{r)vn{r) 



2k, 



-[t]* ■ Vun + V • {unT]*)] = e„w„(r). 



(2.4a) 
(2.4b) 



The wave functions are normalized in two dimensions: / dx J dy 



1. When {un, Vn) is an eigenfunction 



with energy e„, (w* , u* ) is an eigenstate with energy 
as 



-e„. With these eigenfunctions the field operators are expanded 



E 

e„>0 



— ient/h 



(2.5) 



The ground state |0) satisfies 7„(j|0) = 0. 

Gapless chiral edge modes are present at boundaries of a chiral p-wave superconductor. To illustrate this within 
our model, it is sufficient to solve the BdG equation with a simplified gap function with step-function form, t] = 
<d{x)A{l,i£), where <d{x) is the Heaviside step function and e = ±1 is the chirality of the condensate. The self- 
consistent solution of the BdG equation will be presented in the next subsection. The bound state solution to 
Eq. (2^) with energy eigenvalue e = Asin6' (— 7r/2 < 6 < n/2) and with boundary condition m„ = ti„ = at a; = is 
given by 



ue{r) 
ve{r) 



■ exp 



X 

To 



iekpysmO ) sin(A:i?a; cos ( 



gi7r/4 
g-»7r/4 



(2.6) 



where Ly is a linear scale of the superconductor in y-direction. We have used the Andreev approximation valid for 
A ^ /I. The amplitude of the edge state decays in the bulk on the length ^0 = Tivp/A. With the wave number in the 
y direction ky = ekp sin0, the energy dispersion of the chiral edge mode is c{ky) — eAky/kp- 

We now calculate current carried by the chiral edge mode at zero temperature. The current density J is defined as 



I d^rUyiF 



Zim mc 



(2.7) 
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where the fields V'o- are expanded as in (2.5). At this level a finite contribution conies only from the edge mode (p.q). 
Thus, we may restrict the summation over the eigenstates to < < 7r/2. The total charge current running along 
the edge is given by 

/•OO /"OO 

ly^ dx{0\Jy\0) = -— dxlm{0\i;ldyi;i +^jldytl;i\0) 
Jo ^0 

= fdisine) r dxlu,{vedyv;) = (2.8) 

where we have set A ~ 0, ignoring the vector potential induced by the spontaneous current, i.e., diamagnetic screening 
currents. The current ly is spontaneously running along the edge of the chiral p-wave superconductor. If the chemical 
potential could be shifted by a constant as ^ ii + eV in some way, then the spontaneous current would change 
by ee^V/h. From this simple-minded argument it is tempting to conclude a universal value of the Hall conductance. 
This argumentation is invalid, because both the superconducting condensate and the edge states carry current and, 
furthermore, the constant shift of the chemical potential is not realistic to describe a Hall measurement. Indeed a 
careful analysis of the physics of the superconductor surface region is necessary as we will demonstrate below. 



B. Solution for Hall bar geometry 



In this section we study the transverse voltage induced by an externally driven current by solving the BdG equations 
( p.4[ ) self-consistently for a system with Hall bar geometry. We model the Hall bar by a two-dimensional system of 
width Lx in the x direction and length Ly in the y direction. The currents run along the y direction, in which we 
impose the periodic boundary conditions. The superconducting state with the symmetry of the chiral p-wave state 
is parameterized by r; = (Aa;,iAj^), where A^; and Aj, are real functions of r. The calculation is done for zero 
temperature. 

In the Hall bar geometry it is convenient to use the following two basis sets of wave functions, 



V'sin(fc,'r) 



LxLy 



sin kxX 
cos kj^x 



(2.9) 



where ky = 217: / Ly and k^ — mn/Lx {l,m : integer) by setting either the wave function or its first derivative to zero 
at a; = and Lx- We solve the problem twice: (i) with the Dirichlet boundary condition at a; = and Lx, and (ii) with 
the Neumann boundary condition. Afterwards we take the average of both solutions in order to remove unphysically 
rapid changes of the electron density at the surface. Expanding Un{r) and f„(r) as 



Unir) 
Vn{r) 



Un{k) 

(k) 



(2.10) 



we can solve the gap equation numerically by diagonalization. The order parameter is determined self-consistently 
from the gap equation (2^) for T = 



(A^,iAj,) 
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B„>0 



K(r-)Vz;„*(r) - <(r)Vw„(r)] 



(2.11) 



where the coupling constant g is chosen to give A^ 
obtained by solving the gap equation: 



1 



{i-Kkpy 



d^k- 



y/el + {Aok/kpy 



Aq in a bulk chiral p-wave superconductor. That is, g is 
1 + e/fi 



de- 



(2.12) 



where lOc is a cutoff energy and 

1 



h = 



de 



= In 



Ao'(l 



2lJc + AoVa* + 2y^cj2 + Ao2(l+a)c/M) 



-2uj, 



+ AoVa^ + 2^j[VA^il^ZJ^) 



de 
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Ao^ 1 



Ao^i 
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The solutions to the BdG equation (vM determine the electron density and current density: 



(r) = 2 5] K(r)p, 

E„>Q 



Jyir) - - E 

m ^ — ^ 



dy 



dy 



-n^{r)Ayir), 



(2.13a) 
(2.13b) 



where only the y component of the current is non-vanishing in the Hall bar geometry. Note that these densities 
are defined in two dimensions because of the normalization condition imposed on u„ and Vn- The scalar and vector 
potentials obey the Maxwell equations, 



47r 
cd 



Jyir). 



(2.14) 



We have introduced the length d to convert the area densities into the volume densities. Physically d corresponds 
to the spacing between two-dimensional layers in Sr2Ru04. To keep overall charge neutrality, we introduce n^, as 
the density of the uniform positive background charge (jellium model). The externally injected current / fixes the 
boundary conditions for the vector potential: 



dA, 



dx 



dA, 



dx 



2tt 

cd 



-I. 



(2.15) 



The s elf-c onsistent order parameters, scalar and vector potentials can be obtained numerically by solving (2.4), ( 2.11 ) 
and (2.14) iteratively until convergence is reached. In the iteration step we fix the total number of electrons to the 
normal state value by adjusting fi. In the self-consistent solution electric charge is screened on the Thomas-Fermi 
length scale, £ = (?i^(i/4e^m)^/^. Since we have three material-dependent parameters, Ag, fc_F, and d, we have freedom 
to change three dimensionless parameters, kp^o, A/^o, and £/^o- 



C. Discussion of the self-consistent solution 



The solutions of the BdG equation reveal that the relevant physics of the Hall bar indeed happens at the two 
surfaces. First we consider the solution for the case where total current along the y direction is zero. The order 
parameter varies strongly at the surface: A^, is suppressed while Aj^ rises slightly as shown in Fig. 0(a). This behavior 
is connected with the reflection properties of Cooper pairs at the surface. Resulting interference effects are destructive 
for Ax, since the order parameter is odd under reflection at a surface normal to the x direction. Note that this order 
parameter variation is specific for a specular reflection at the surface and would look different for the case of diffuse 
scattering. We do not, however, consider this aspect here further. 

At the surface the chiral edge states appear with a linear dispersion around the Fermi energy. These are Andreev 
bound states as a direct consequence of the chirality of the superconducting state. The two branches seen in Fig. |l|(b) 
belong each to one of the two edges of the bar. These chiral edge states generate spontaneous currents at the surface 
flowing along n x z (n: surface normal vector). The currents on the two edges of the Hall bar run in opposite 
directions, thereby no net current flowing in the Hall bar. The edge currents generate a magnetic field, which is 
screened by counter currents in the interior of the superconductor [Fig. |l|(c)]. The length scale of the surface current 
is the coherence length while the screening currents spread over the London penetration depth A. The surface 
magnetic fields generate a finite magnetization whose sign depends on the sign of Cooper pair angular momentum 
and the sign of the charge, i.e., electron- like or hole-like Fermi surface. 

Turning now to the question of the scalar potential and the charge distribution in the superconductor, we find again 
that all interesting features show up only in the surface region. There is a finite excess charge at the surface, which 
is screened due to Thomas-Fermi screening. As a result the charge density forms a dipole layer and is overall charge 
neutral. This is the constraint that we have imposed by fixing the external electrical field to zero (note that in our 
Hall bar geometry a finite charge or equivalently a finite external electric field would correspond to an infinite field 
energy, since the capacitance is zero). The dipole layer induces a local electric field and causes a shift of the scalar 
potential relative to its value in the bulk of the superconductor which we choose to be zero. The charge distributions 
and potential at the surface are the same on both sides of the bar. 

Under the assumption that the Hall bar is symmetric about x = there is no potential difference between 

the two sides. We would like to mention that scalar potential variations close to the surface are not unique to chiral 
superconducting states, but occur in any superconductor whose order parameter is influenced by surface scattering. 
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If we introduce a net current driven from an external source, this external current will be distributed equally to 
both surfaces [see Fig. §(a)]. This affects the surface states differently for the two sides, because on one side the 
external current flows with the spontaneous current, on the other against it. Furthermore, the scalar potential is not 
equal anymore at the two sides because the charge dipole layers are modified differently, as we can see in Fig. ^(b). 
This transverse voltage difference depends on the orientation of the external current and appears in the absence of 
an external magnetic field, thereby referred to as the SHE. 

Our calculation clearly shows a linear relation between the source current / and the transverse voltage Vh as 
expected for the Hall effect. Deviations occur only when the current approaches the critical value where the order 
parameter starts to be strongly affected by the current. In Fig. || we show the k, i and kp^o dependences of the 
Hall resistance, Rh = Vh/I- What is immediately obvious is that the Hall resistance is strongly suppressed from the 
quantum unit of resistance Rq = /i/2e^. There is also a strong dependence on the material dependent parameters, 
indicating that the behavior is non-universal. Unfortunately, in the numerical BdG scheme we are limited in the 
choice of i, ^o, A and kp^, because large difference in their magnitudes demands large-scale computation. In the 
next section we study the SHE using the extended GL theory which will allow us to calculate Rh analytically for 
temperatures close to Tc- We can already here confirm that the quantitative comparison between the two methods 
work very well. It will also become clear that there are a few contributions to the transverse voltage. 



III. GINZBURG-LANDAU FORMULATION 



We formulate an extended Ginzburg-Landau theory based on symmetry arguments which includes the scalar and 
vector potentials in a general form. This allows us to analyze the anomalous coupling between charge and magnetic 
degrees of freedom in chiral superconductors. 



A. Ginzburg-Landau free energy 

The paring symmetry of the chiral p- wave superconducting state is characterized by d(k) = zA(kx±iky)/kp, which 
is a combination of the two degenerate p-wave components with p^- and py-symmetry. This degeneracy is not lifted 
when we introduce a tetragonal crystal field, although the details of the fc-dependence of d may change. The two-fold 
degeneracy requires that we introduce two complex order parameter components, r] = {rix,ily), such that the d- vector 
becomes d{k) = z{ri ■ k)/kF- The free energy has to be a scalar under the transformations of the symmetry group 



g = DihxTx U{1), 



(3.1) 



where the tetragonal point group D^h includes the simultaneous transformation of orbital and spin degree of freedom 
as a consequence of spin-orbit coupling, T denotes the time-reversai-symmetry and U{1) the gauge symmetry. The 
GL free energy expansion for t] with the symmetry Q is well known:E3 F = J <PrJ-, where 

+Ki {\D,7^x\^ + \DyT]y\^) + K2 {\D,Vy\^ + \Dyri,\^) 

+ K3 [{D^n^YiDyl^y) + C.C.] + [(Dxl]y)*{DyJ]x) + C.C.] 

(V X A)2 



(3.2) 



The coefficients a, bi and Ki are non-universal real numbers that depend on the details of the material. The coefficient 
a is negative below Tc (a cx T — Tc). For the choice < ^2 < 4fei + 63 and 63 < 62, we find the homogeneous phase 
r] = 770(1, ±i) with ?7q(T) = |a|/(46i — 62 + b^). The gradient terms are explicitly gauge invariant by the definition 
D — V + i{2e/hc)A. Equation (|3.2| ) is the standard free energy density used to study the response to the magnetic 
field with the Coulomb gauge V • A = 0. 

As is well known,cil the intrinsic orbital angular momentum in a chiral p-wave state is related to the difference of 
the two terms proportional to K3 and K4: 



{DxVx)*Dyriy - (DxTiy)* DyTjx + c.c. 



V X {r]*Dr,y - r,lDr]x) 



r];r,y)V X A 



(3.3) 
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where it is understood that the z component of a three-dimensional vector is kept in the right hand side. The integral 
of the first term gives a surface term which we may discard, while the second term represents the Zeeman energy of 
a magnetic moment coming from the intrinsic angular moment (oc iry x r/*). 

For our purpose it is essential to include additional terms that are coupled to either the scalar potential or 
the electric field E — — V^o- We consider the stationary situation and fix the gauge so that Aq is zero in the 
undisturbed homogeneous superconductor. It is more convenient to use the Lagrangian formulation, in which the 
Maxwell equations simply follow from variation with respect to A^. The Lagrangian density involving the scalar 
potential is given by 

+K3Ao[{Dxr]x)* {Dyr]y) + {D^T]x){DyT]y)*] 

+K4Ao[{D^'ny)*{Dy'n^) + {D^'ny){DyVxT] + KsAoHD^r^^f + \D,T]yf) 
+Ci[E^T^l{D.^ri^) + EyT^liDyTiy) + C.C.] + C2[E^vl{DxVy) + Eyii^DyVx) + c.c] 
+Cz[E^T^l{Dyriy) + EyT^l{D.^ri^) + c.c] + C4[E^r]l{DyT]^) + Eyr]*{D^T]y) + c.c] 

+C,[E^V;{D^V.) + E.v;{D.Vy) + c.c] - - ^ (3.4) 

with Ki,Ci and £ being real numbers. The coefficients will be derived for the weak-coupling limit in Sec. V. In this 
form it is easy to verify that each term is individually invariant. The term Aq/8tt£'^ describes the screening of the 
electric charge in the metallic and superconducting state, where £ is the Thomas-Fermi screening length. The choice of 
this form fixes Aq to zero in the bulk of the superconductors, which corresponds to the chemical potential as required 
by our choice of gauge. We emphasize that Ce + T plays a role of the Lagrangian density for A^. 

We notice that the Ki terms are closely related with the Ki terms. On the other hand, the Ci terms have no relatives 



in Eq. (3.2). The difference of the C3 and C4 terms contains a contribution i{rjx'ri* — r]y'q*^){AxdyAt^ — AydxA^)). It 
describes a coupling between the scalar and the vector potentials and is similar to the Chern-Simons (CS) term. Since 
the fields are static, however, it is not exactly the same as the CS term. Thus we shall call it a CS-like term. It can also 
be interpreted as representing a coupling between the intrinsic magnetic moment (oc irj x rj*) and the magnetic field 
in the presence of an electric field. In other words, the CS-like term describes the reaction of the intrinsic magnetic 
moment to a change in the Cooper pair density.Ej The CS-like term proportional to C3 — C4 will play an essential 
role in the SHE besides others. 



B. Equations for the electromagnetism 

We derive equations describing the electromagnetic properties of the superconductor from variation of J (Pr(Ce +^) 
with respect to Aq and A. The equation for the scalar potential has the form, 

-V^Ao + ^=4n[p-V-{P + n)], (3.5) 



where 



p = Ki{\Dxi^xf + iDyT^yf) + K2{\DxVyf + iDyV.f) + MiD^ij^)* {Dyi^y) + {Dxi^x){Dyinyr] 

+K4iD,r]yyiDyrj,) + {D,r^y){Dy7^,y] + K^HD^t^^I'' + \D,r]yf), (3.6a) 

'd,iCi\rj,\^ + C2\riy?) + i(C3 + Ci)dy{v;Vy + V.V*y) ' 

dy{Ci\ijy\^ + C2\v.?) + U'^3 + C4)dxi'n*Vy + ijxT^;) |, (3.6b) 



1 ^ _ /' ^y^y'^^ + ^v^y"^^ ^ V*xDyVy - VxD*jj*^ 
2*' 



n = - (C3 - C4) v*xD,Vy + VxDlv; - V*yDxVx - VyD*V*x ■ (3.6c) 







The two terms on the left hand side of Eq. (1_5) describe the screening of the electric field. The quantity p is the 
charge density induced by variations of the order parameter. The vector density P represents the electric polarization 
caused by inhomogeneities of the superconducting condensate. Both p and P appear as a consequence of the fact 
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that the local change in the condensate, both in the order parameter and in the supercurrent, leads to a redistribution 
of the electric charge. This is most easily seen in p, which is related to the gradient term of the free energy. Note 
that terms of this kind are also present in conventional superconductors. On the other hand, the vector density 11 is 
anomalous and characteristic of chiral superconductors. 

For the sake of simplicity we assume that the relative phase of the order parameter component is fixed, i.e., 
Vx = IVxl exp(i(/)) and rjy = *£|?7al exp(j(/)) with e = ±1, where e is the chirality of the p-wave order parameter. This 
condition is satisfied in the situations we will study below. In this case we have i(r]*r]x — 'ri*r]y) ~ 2e\rix\\riy\. The 
vector n can be written as 



2e 

he 



Its divergence is equivalent to a source charge, 

4e 



Pn 



V n 



^^e\7^x\\Vy\B,{C3 - C4) - 2e{C3 - Ci)iV\rjx\\Vy\) x I V0 + 



2e 



(3.7) 



(3.8) 



The first term indicates that the magnetic field B induces the electric charge whose sign depends on the chirality e. 
The second term is non-vanishing only when the modulus of the order parameter has spatial variation. Terms similar 
to this also exist in p. 

The modified London equation is obtained by variation of J cPr{Ce + T) with respect to A: 



V^A+^ (j + J + y) 



0, 



where we find three current contributions. The first two, 



4e 



and 



J 



'Im {KirilDxrix + K2r]*Dxriy + K'^rjlDyTiy + Kiri*Dyrix 
J = - ^ ( Im {Kii]*yDyrjy + K2ri*Dyr]x + KsrjyDxTjx + Ki-qlDxily) 

Im {K5r]*D^rix + K^-qlD^riy) 

^Im {KiTjlDx'qx + K2i]*Dxriy + Kj,rilDyriy + Ki-q^Dyij^ \ 
Im {Kir]*yDyriy + K2rixDyrix + K3ri*Dx'nx + K4J]*Dxr]y 
Im (Kr,r]*D^r]x + Kr^rflD^'qy 



(3.9) 



(3.10a) 



(3.10b) 



are the supercurrents including the screening currents. The second one contributes only if the scalar potential deviates 
from its bulk value Aq — Q. The last term 



2e 

Y = i^(C3 - C4){VxV*y - VxVy)E X z 



(3.10c) 



is the anomalous contribution, where the electric field acts as a source of supercurrent. Note, however, that Y 
does not cause dissipation because E ■ Y = 0. Both anomalous components 11 and Y are proportional to the 
chirality £ and originate from the CS-like term. We would like to emphasize here that the presence of the anomalous 
current contribution does not affect the standard flux quantization. The current Y quickly drops to zero inside the 
superconducting material because the electrical fieldpis strongly screened. Thus even a hole containing a net charge 
would not violate standard laws of flux quantization .EJ 



IV. TRANSVERSE VOLTAGE IN THE GL DESCRIPTION 



Now we shall study the SHE in a Hall bar geometry on the level of the extended GL theory. Since we concluded in 
the BdG study that the relevant physics lies in the behavior of the superconductor close to the surface whose influence 
is exponentially small in the interior, we concentrate only on one edge and consider the half-space x > 0. The physics 
of a Hall bar of a given width follows in a simple way from the results obtained for a single edge. 
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A. The surface state 



The boundary conditions for the order parameter are chosen assuming specular surface scattering. For the given 
geometry they read 

?7xU=o = and Dj^ijyl^^o ^ 0, (4.1) 

implying that there is no current running normal through the surface.0 We need to consider only the x-component 
of the electric field and the z-component of the magnetic field, both of which are continuous at the surface and are 
functions of x. It can be easily shown that these boundary conditions are compatible with the gauge invariance of 
J-e- The self- inductance L and the capacitance C of our system are taken to be infinity and zero, respectively, such 
that net current and total charge in the system should vanish, unless they are imposed by external sources. This is 
important for the choice of the boundary conditions for the gauge fields. 

We present an approximate solution to the GL equation that captures the essential aspects of the problem. We 
assume that the system is in a chiral p-wave state of a single domain with the chirality e. First we solve the GL 
equations to determine the spatial dependence of the order parameter r]. To this end we ignore the gauge fields A^^ 
in the GL equation. This is allowed because we are looking for a solution in lowest order in t = {Tc — T)/Tc, a small 
parameter in the theory. The GL equations for rjx ~ \rix\ and rjy = ie\T^y \ are 

a|?7,| + 26i|77,|3 + (26i - 63 + h)\Tj,\\i^y\^ , (4.2a) 

a\riy\ + 26i|?7y|3 + (26i - 62 + h)\n,\''\i^y\. (4.2b) 

To simplify the analysis we consider a special situation where the coefficients hi in the GL free energy satisfy the 
relation 



K2 



\Vx\ 



<P\Vy\ 



2bi = &2 - b3, 



(4.3) 



which, however, is n ot sat isfied for a cylindrical Fermi surface in the weak-coupling limit (see Sec. V). Under this 
condition Eqs. ( 4.2a| ) and ( 4.2b ) have the solution. 



7/0 tanh 



Vy = iSVo 



(4.4) 



with ?7q = |a|/26i and the c ohe rence length = 2Ki/\a\. Note that rj^ vanishes linearly at a; = whereas r]y stays 
constant. If the condition (4.3) is not satisfied, rjy deviates slightly at a; ^ ^ from the bulk value rjo, but generally 
does not vanish at x = 0, as we have seen in the BdG study. In general the order parameter shows the following 
feature: the component of r] normal to the surface vanishes linearly while the perpendicular component is only weakly 
affected. As the detailed spatial dependenc e of r) is not important for our semi-quantitative discussion, we will use 
an approximate form for rj^, instead of Eq. (4.4), which allows us to proceed with analytic calculations more easily. 



Vx 



(4.5) 



We would like to mention that we have solved the coupled GL equations for more general cases by numerical means 
to confirm that our approximations work very well both qualitatively and quantitatively. 

Having determined the profile of the order parameter, we can calculate the distribution of the charge and the 
spontaneous current in the equilibrium state from Eqs. ( |3.6D and (3.10). The equations for the scalar and vector 
potential have the following form: 



dx^ £2 



d^A. 



y 



A 



y 

dx"^ ^1 V J dx 6^3(2:) ^A^4 

where we have introduced parameters h {i — 1 ^ 4) that have dimension of length: 



dAy 
dx 

dAp 
dx 



eA, 



2AoAy 



6/3(2;) 



2Cie- 



(4C1 + ifi) 



-2x/i 



(4.6) 
(4.7) 
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/ 


= (^^^^ 


1 

7 




1 




K 








1 




3 


h 


he 





The London penetration depth A is given by 



1 

A2 



327re^ 



iKi+K2). 



(4.8a) 
(4.8b) 
(4.8c) 
(4.8d) 

(4.9) 



In deriving Eq. ( [4.7D we have replaced 



1 

A2 



by ignoring the spatial dependence. This approximation weakly affects a numerical factor in the final expression 
of the Hall voltage. 

For temperatures close to the onset of superconductivity, r — {Tc — T)/Tc is a small parameter which allows for a 
controlled approximation. From the standard temperature dependence of 770, ^^'^ ^^ it is clear that 



(4.10) 



We first consider the case of vanishing external fields, i.e., E = B = in the vacuum (x < 0). This means that the 
net charge and current in the superconductor are zero. In this case the scalar potential Aq""* and the vector potential 



A'y' obey the boundary conditions dA^Q'/dx = dA^y'/dx = at a; = 0. From Eq. ( [l.lOD we see that A^^' = ©(t^) 
and A^y^ = 0{t^^'^) since £ = 0{t^) and ^ ^, A. In lowest order in r we find 



4°)(x) 



4°^(x) 



dAf\x) 



(4.11) 



1 



^(3Ci+i^i)e- 



dx 

2Cie"^/« - (4C1 
where k = X/^. This gives the total magnetic fiux per length, 

<i>(°) = / dxB^ix) = -4°)(0) = -ef ^ 

Jo liK+1 



.4°)(0)e 



6/3(2:) 



(4.12) 



(4.13) 



Note that the vector potential A'"y^ is proportional to e whereas the scalar potential is independent of the chirality. 
For both the scalar and the vector potential the London penetration depth constitutes the longest length scale of 
variation, because they are coupled together. The charge density —{l/4:Tr)d'^Al^''/dx'^ has dipolar form as found in 
the BdG calculation. The total current density is 



.(0) c d'^A^y\x) ce 
■'y ~~4^ d^ ~ ^AttUX^ 



-^(«e-«-.-A) 



(4J4) 



where the first term is related to the spontaneous current due to the chiral Andreev bound states, and the second 
term is the response of the superconductor, i.e., the screening current. Thus the net current is zero in the absence of 
an external magnetic field. 

Results obtained for the semi-infinite geometry can be easily carried over to a Hall bar that extends from x — 
to X = L^. When the width of the Hall bar is much larger than A, the two surfaces are basically disconnected 
electromagnetically. The gauge fields for the Hall bar geometry are then obtained by combining the contributions 
from the two edges. The scalar potential is A^^ (x) + A^^ (L^ — x) , while the vector potential is Ay"' (x) — A'"y^ {L^ — x) . 
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B. Response to external fields 



We now consider two cases where a weak external field is applied to the semi-infinite system. Weak perturbations 



introduce corrections to the gauge fields, Aq — A'j^^ + SAq and Ay = Ay'' + SAy. We would like to obtain 5 A 
linear response to the external perturbation. We can, therefore, linearize the Maxwell equations in 5Aa: 



(0) 



efSAo 5Ao e 



(f6Ay 



5Ay 

A2 



dx 
dx 



-(- 

e 



6^3(2;) 



2A 



(0) 



6^3(0; 

[A'-^hAy 



■SAy = 0, 



(4.15) 
(4.16) 



(1) In the first case a weak external magnetic field Bz is present in the vacuum, which corresponds to a finite net 
current running in the y direction. The boundary conditions are dSA^/dx = and dSAy/dx = B sd x = 0. In leading 
orders in r we obtain 



SA<y'\x) 

S4'\x) 



-eBf 



(4.17) 



1 1 

K 



-X j \—x I ^ 



-x/X~x/i 



-x/l 



.-x/i 



h{0) 



(«:-l) 



(4.18) 



We note that the sign of the scalar potential in linear response to B depends on the chirality e. 

This result can be used to determine the transverse voltage in the Hall bar. We assume that the induced current 
fiows symmetrically on the two edges. Formally this situation is realized by applying the field Bz[Q{—x) ~Q{x ~ L^)]- 



The change of the scalar potential is SAI^\x) ~ SA^^'iL^ — x), leading to a finite transverse voltage across the Hall 

bar, Vh = Ao(Lx) — Aq{0) = —25A'^\o). Now we can make connection to the SHE in a two-dimensional system or in 
a layered system like Sr2Ru04, where the system consists of layers (xy-planes) separat ed by a distance d. The total 
current per layer / is related to the external magnetic field by / = cdB/2TT. From Eq. ( 4.18 ) we find 



1(1) 



VH=eI 



47r^2 
cd 



2A 



;3(0)/4 K + 1 



(4.19) 



where we have kept only the leading terms, which are proportional to t. The parameters entering this expression will 
be derived from a microscopic model in the weak-coupling limit in the next section. 

(2) We now turn to the second case where a weak external electric field Ex is applied to thej-superconductor, 
inducing a finite surface charge. A similar problem was recently discussed by Goryo and Ishikawa.li3 The bou ndary 
conditions here are ddA^ j dx 

and dTTel), 



-E and dSAy/dx = at x = 0. In lowest orders in r we find the solution to Eqs. (4.15) 



SA^^Hx) ^ -eEfe--/' 



2£ 



2£ 



+eEte 



-x/X 



1 1 



2A K 



(4.20) 



(4.21) 



23(0)^4 K+1. 

The external electric field changes the local configuration of electric current and magnetic field in the chiral super- 

(2) 

conductor. The induced vector potential 6 Ay depends on the chirality e. The total change in the magnetic flux (per 



unit length along the y direction at the surface) is 

(5$ = -SA'i^^O) = eEe 



1 1 



2A 



^3(0)^4 K + 1 



(4.22) 



Note the similarity between Eqs. (4.19) and (4.22) indicating that the two phenomena have indeed common origin. 



11 



V. THE WEAK-COUPLING COEFFICIENTS OF THE GL THEORY 



In this section we calculate the coefficients of the GL free energy in the weak-coupling limit using a model with 
layers of two-dimensional electron gas, where electrons are confined in each layer. Furthermore we ignore any spatial 
variation of electromagnetic fields in the direction perpendicular to the layers. 



T he m ean- field Hamiltonian appropriate for the discussion of the chiral p-wave state (1.1) is already given in 
Eq. {2A). The static electromagnetic fields are governed by 



Hem = ^ [(V X Af - (V^o)'] 



(5.1) 



where d is the distance between the two-dimensional layers. The sign of the second term was chosen negative so that 
the Maxwell equations can be obtained by taking the functional derivative of J d^r {TCmf + '^em) with respect to A^, 
which should be viewed as the Lagrangian in imaginary time. 

We integrate out the electron fields ipa and ipa to obtain the effective functional i^off for t] and A^: 



Yl Vi^^V^P^ exp - - y dt j d^r {hi^^dti^a + Umf + Hem) 



(5.2) 



where the electron field operators '0cr(r) and '(pti'"') in Wmf (2.1) are replaced by Grassman fields ipaifjit) and 
ipaifjit), respectively. The GL equations are then obtained by taking the functional derivative of Fcs with respect 
to T] and Afj^. We calculate -Fcff in powers of iji and up to the order rjf, DirjjDkrji, AoDirjjDkrji and EirjjDkrji in 
the weak-coupling limit, jr/j ^ h^kp/2m. The calculation is tedious but straightforward, and only the final result is 
presented below. The functional Fos so obtained has the form 



d^r {T -f Ce 



(5.3) 



The free energy part JF, which is defined in Eq. (|3.2| ) as the free-energy density in three dimensions, has the standard 
coefficients, 



3 



iV(0), 



63 _ 7C(3)iV(0) 



2 128(7rfcBTc)2 ' 

7C(3)A^(0) / hvF 



128 



(5.4a) 
(5.4b) 

(5.4c) 



where depends on the coupling constant g in the usual exponential form ksTc — LUcexp[—2/ gN{0)] with lUc being 
the cutoff energy scale, and N(0) = m/{2iih^d) is the density of states (per spin) at the Fermi level. The results (5.4) 
are valid in lowest order in |j7|//i, where we find = K^. If we assume that the density of states had a weak energy 
dependence with energy derivative iV'(O) ~ N{0)/fi — / {irK^kpd) at the Fermi surface, then there would be a tiny 
difference between K3 and K4, yielding the contribution to the free energy density. 



8cm 



2e*^WD 
irTr 



[iT]{r) X 77* (r)] 



2^2 



B, 



(5.5) 



where C = 0.5772 ... is the Euler's constant, uq = kp/2TTd is the electron density, and /x is the magnetic moment 
(per electron) of a Gooper pair in the chiral p-wave state. 



(5.6) 



Here iib = eh/2mc is the Bohr magneton. Enuatinn (5.5) is the Zeeman energy for the intrinsic magnetic moment 
of the chiral p-wave state. As is well known,t3L3l£3 it is diminished by the factor (770//^)^ which indicates the degree 
of particle-hole asymmetry at the Fermi level. Since this contribution is very small, we can ignore it in the following 
analysis. 



From the coefficients in Eq. (5.4) we immediately obtain the order parameter. 
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ryo(T) = nksTc 



I 8r 
7C(3)' 



the coherence length, 



and the London penetration depth, 



A(r) 



he 



2Ki hvF /21C(3) 




irksTc V 32r 



2er]o V 87r(i^i + K2) 




8tttN{0) 



The second contribution to the free energy, Ce, has the foUowing coefficients: 

1 



87r€2 



3 1287rd(7rfcBTc)2 ' 



C3 



C4 



7C(3)e 



2567rd(7rfcBrc)2 



(5.7) 
(5.8) 

(5.9) 

(5.10a) 
(5.10b) 
(5.10c) 



Notice that Ki = edKi/ dfi. T his rel ation can be easily understood, if we regard eA^ as a spatial variation of /i. In the 
course of deriving Eqs. ( ^.10b| ) and ( ^.lOcD we have naturally assumed that the momentum k in d{k) = z{r] ■ k)/kF 
is close to the Fermi surface. The CS-like term describing the reaction of the intrinsic magnetic moment to a change 
in the superfluid denaiUi; does not have the reduction factor found in Eq. (5.5), in agreement with the argument by 

Volovik and Mine evBjg 

From Eqs. (5.4) and ( 5.10| ) we obtain the parameters appearing in Eqs. (4. IS) and (4.22): 

hcd 

(5.11a) 



' 2 " 2e2r' 



A 



^3(0)^4 



ihcd 



(5.11b) 



Now we may express the results for the two cases discussed in the last section in terms of the microscopic parameters. 
The Hall voltage induced by an external current is found to be 



Vh = el 



16(efcFA)2 



2n 



3k 



(5.12) 



in the leading order in r. Here we have made use of the relation (tpX)'^ — mc^d/Ae'^T. We find a strong reduction 
compared with the quantum unit of resistance i?o — h/2e'^. The factor l/{kp\)^ can be also written as IGnnsx/nQ, 
where Ug is the su perfluid density and x = A'b^(O) = 3xorb corresponding to the orbital susceptibility. We can thus 
rewrite Eq. (5.12) as 



Vh 



sl^{2n- 



3k 

K+1 



(5.13) 



where 7 = tIs/uq — 2r is the ratio of the superfluid density to the electron density. Obviously the Hall resistance 
strongly depends on material dependent parameters and is also rather severely reduced from the universal value h/e'^. 
In Sec. VI we will analyze the quantitative aspect in more detail. 

We turn now to the reciprocal case (2). Here the effect is more subtle as the response to an external electric 
field constitutes a change of the field distribution in the vicinity of the surface. Since there is already spontaneous 
magnetization generated by the surface currents, we estimate the relative change of the total magnetic flux: 



(5$ 5Ai^\0) 



$(0) 



4°'(o) 



(5.14) 



The electric field E should be smaller than A/e£ in order to avoid nonlinear effects that arise from the field effect on 
the superconducting order parameter at the surface. Therefore the induced change of flux (5$ would only be a very 
small fraction of the spontaneous magnetic flux 
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VI. DISCUSSION BASED ON HYDRODYNAMIC EQUATIONS 



A. The hydrodynamic equations 



In the previous sections we have considered the SHE using the self-consistent solution of the BdG equation at 
T — and the extended GL theory near Tc. In this section we introduce a phenomenological description based 
on the stationary hydrodynamic equations. This approach can provide an interpolation between the two limits and 
allow us to have a simple intuitive understanding of the physics involved. Our starting point is the phenomenological 
Lagrangian 



OTT 



(6.1) 



which describes electromagnetic properties of a superconductor. In the chiral p-wave state the Lagrangian density / 
may be written as 



/ 



(6.2) 



where Vg — {e/mcjA and Vq is the velocity of the supercurrent generated by the spatial dependence of the order 
parameter: 



Vox 

VOy 



8m?7Q 



Im [i-qldx-qx + VydxVy) 
Im {ir^ldyT^y + vldyi^^) 



f Im {vldyi^y + ijydyrjx) 
vim {'n*d^,r]y + r]*d^r]^) 



(6.3) 



The first term corresponds to the ordinary supercurrent due to phase gradient, while the second term is connected 
with the spontaneous current due to texture of the order parameter. The latter is equivalent t o th e surface current 
of the chiral edge states and will be important for our discussion. The partial derivative in Eq. (O) acts on rig as 



971s dus dnn „ 



(6.4) 



p/2TTd. We know the two 



where 7 = rig/ no is the ratio of the superfluid density Ug to the electron density no = 
limiting values of 7: 7 = 1 at T = and 7 = 2r near T^.. 

The Lagrangian density / can be deduced from the GL Lagrangian in the following way. In the chira l p-wave state 
we may assume without loss of gene rality that r]*r]y is imaginary. With the weak-coupling result ( ^.4c ), it is easy to 
confirm that the Ki terms in (3.2) generate UgrnVg -Vq, in addition to ns{m/2)v'^. The latter term can be obtained 
with the approximation \r]x\'^ = Ivyl'^ = Vo- then natural to complete the square to make m(vs + Vq)'^ /2, which 
is nothing but the kinetic energy of the superfiuid in the presence of the spontaneous fiow with the velocity vq- From 
the relation Ki = edKi/d^, one can also s ee tha t the Ki terms lead to the kinetic energy contribution to eAodf /dfi. 
Furthermore, we find from Eqs. (3.4) and (5.10c) that the Ci terms have a contribution 



8e f 

-«C'i^ / d^r{r]lriy - 7j*ri.^){AydxAo - Ar^dyAo) 
8e f 

= ^C*!— / (frAoV ■ [a X z{T]lr]y — 77*771;)] -I- surface term 



(6.5) 



in the chiral p-wave state where 77*77j^ is imaginary. This gives the remaining term, 2^eAQN (0)11^^ ■ (/x x A). With 
the help of the identity V ■ {fj, x A) = A ■ {V x fi) — fi ■ {V x A) , we find that V ■ {fi x A) represents the coupling 
of the magnetic current V x /j to the vector potential as well as the Zeeman energy of the magnetic moment fj,. It is 
impo rtan t to realize that the full magnetic moment of a chiral Cooper pair appears here, without reduction found in 
Eq. (5.5). The scalar potential induces a change in the nuaiber of Cooper pairs, which is necessarily accompanied by 
the change of the full magnetic moment per Cooper pairnJ 

We now take the variation of F with respect to A to obtain the extended London equation. 



V'A 



4Tr 



0, 



(6.6) 



with the current density 
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j = -UsCVq - 2ejN{0)[eAo{vs + vq) + cE x /x]. 



(6.7) 



The variation of F with respect to Aq yields 



£^ no 



(6.8) 



Since the Thomas-Fermi screening length £ is much shorter than ^ and A, we may ignore V^Aq to obtain immediately 



An 



eno 



V -ipx A) 



(6.9) 



which we may call a generalized Bernoulli equation. The first term represents a Bernoulli force coming from the kinetic 
energy of superfluid. In conventional superconductors without broken time-reversal symmetry, the spontaneous current 
is absent. In such a case a supsrpurrent / injected from an external current source can generate a transverse potential 
difference proportional to /-^Ea^E^I On the other hand, in the case of our interest where a spontaneous current flows 
along a bounda ry { vq ^ 0), an external current can induce a transverse potential difference proportional to /. The 
second term in (S.9) is characteristic of the chiralp-wave state where Cooper pairs have their magnetic moment. Both 
terms are important in the SHE and in its reciprocal effect. 



B. Spontaneous Hall effect and its reciprocal effect 



Let us study the spontaneous Hall effect using the hydrodynamic equations. As in Sec. IV, we consider a chiral 
p-wave superconductor with a boundary at cc = only. The superconductor occupies the positive x region and 
the system is translationally invariant in the y and z directions. We suppose that the spatial profile of the order 
parameter T]{r) = (r]x{x),'r]y{x)'^ is already determined self-consistently in the equilibrium state. In particular, as 
we have observed in the BdG calculation in Sec. II, the order parameter satisfies ij^ = drjy/dx = at a; = under 
the assumption of the specular reflection at the surface. We will calculate a linear response to week electromagnetic 
perturbation. 

The local variation of 77 near the boundary yields a spontaneous current due to the order parameter texture, 
jt = -Usevt-, where 

r^d^''^^^ 



-t = ^lm{,:^+,;^]y. (6.10) 



This current determines the vector potential A — (0,^j^,0) obeying the London equation 

l \ . 47r 



Jx^ Xy^--T^- ^'-''^ 
whose solution under the boundary condition dxAy = at a; is 

A(0)(a;) ^^J^ (e-l---'l/^ + e-l-+-'l/A^ jt{x')dx'. (6.12) 

In particular, its boundary value is 

Having determined r) and A, we are now ready to obtain the scalar potential ^0 from the generalized Bernoulli 
equation. 

Suppose that a small current / is injected from an external source to the superconductor, yielding a small change 
in the vector potential: A — A^") + SA, where SA{0) — —I{2nX/cd)y. Accordingly, the scalar potential acquires a 
small change, whose boundary value is 



Mo(0) = ^SA ■ ( ^A(o) + -vt + V X 

eno \TOc^ c / 



(6.14) 

x=0 
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where we have used the fact that /i = at x = 0. The intrinsic magnetic moment fi defined in Eq. (5.6) gives 



V X u - -^Im (^^77 +n*^\v 

The spontaneous H all vo lt age in the H all b ar geometry can be related to SAq by the relation Vh 
Sec. IV. From Eqs. ( |6.10|) , (|6.12| ), and (|6.15|) we obtain 



Vh 2-KhnsX 



With the approximate form (4.5) of the order parameter used in Sec. IV, Eq. (6.16) becomes 



Vh_ 
I 



27r/i 



K+ 1 



(6.15) 

-25Aq{Q) as in 

(6.16) 
(6.17) 



which compares well with (5.13) 



The quantitative comparison with the self-consistent solution of the BdG equations encounters some drawback from 
the fact that for the sake of numerical accuracy the cutoff energy uuc had to be taken comparable to the Fermi energy. 
Therefore, there are larger strong-coupling corrections not included in our phenomenological analysis. Nevertheless, 
we find a reasonably good agreement between the phenomenological and the BdG estimates. The result from the 
numerical BdG calculation {kp^a — 16 and ^/^o = 0.1) is 
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while Eq. (6.17) with the same parameters (apart from lo^ leads to 
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K = 1, 

n = 2. 



(6.18) 



(6.19) 



The discrepancy is not only a result of weak-coupling versus strong-coupling approach, but we would like to remind 
that we have also used an approximate description of the order parameter texture at the surface. 

Obviously the resistance obtained here is considerably smaller than the universal unit hje^ . We would like to build 
now a connection between the SHE and the ordinary Hall effect. We may look for the intrinsic magnetic field which 
causes the Hall response to the externally induced current /. Comparing Eq. ( |6.17 ) with the standard expression of 
the Hall effect. 
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we obtain the effective magnetic field 
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(6.20) 



(6.21) 



The effectiv e fie ld corresponds to the density of magnetic moments of Cooper pairs which is not the reduced magnetic 
moment in (5.5), but rather the full moment which is associated with the change of the Cooper pair density. Note 
that the factor 7rK[l -t- 1/(1 -I- k)] is due to inhomogeneous field and current distribution. 

Next we consider the reciprocal effect. In the presence of a weak external electric field SE = Ex at a; < 0, the 
scalar potential in the superconductor receives a small perturbation, SAq{x) = Ele~^/^. This yields a small change 
in the current density. 
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(6.22) 



which determines the change in the vector potential bA. Solving the London equation, we find 
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(6.23) 



Since the Thomas-Fermi screening length £ is much shorter than ^ and A, we may set vty{x) = vty(Q), A'"y\x) = A^y\0), 
and iJLz(x) = x^'^(Q) in 5jy(x) in the integrand. With this approximation we obtain 
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(6.24) 



The ma gnitude of the ratio dAy{0)/A'j^\o) may be estimated from Eqs. ( |6.13| ) and ( |6.24D for the approximate form 
(4.5), yielding 



eE£^ 



A^\0) AiA 



(« + 2), 



(6.25) 



in good agreement with the GL analysis. Unfortunately, in this case the comparison with the BdG result does not 
agree well, which we attribute to the fact that £ is comparable to the other length scales in the numerical calculation. 
Thus, there are strong correction to the above result in addition to the strong-coupling corrections. 



C. Experimental probe 

Although the phenomenon we discuss in this paper can be compared to the standard Hall effect, it has actually 
some quite distinct aspects. Flowing currents are an equilibrium property of the superconductor. It is, therefore, 
impossible to measure the transverse voltage by means of a standard voltmeters, using direct contacts to the surfaces. 
This problem was realized alreadv-4nore than 30 years ago, when the Bernoulli response to a current was investigated 
in conventional superconductors J23 In this case a potential difference between the surface and the interior of the 
superconductor is expected, which is proportional to the square of the running current: 

^ - (6.26) 

The potential difference is actually independent of temperature as has been observedEl for Pb. The method of 
measurement was based on a thin-film capacitor which picks up the voltage signal caused by an ac-current on the 
superconductor surface. The same kind of capacitor technique could also be used to detect the spontaneous Hall 
effect. In an ac-measurement the above Bernoulli force (oc I^) yields the second harmonic of the applied ac-current, 
while the signal corresponding to the SHE (cx /) contributes to the first harmonic. Hence we can distinguish the SHE 
from the standard Bernoulli effect. We can estimate the magnitude of the Hall response 

where /tot is the total current running through a three-dimensional sample consisting of Mayer layers. Using typical 
values for ordinary metals (for example, kpX « 200 at T = 0), we obtain Rh ~ 0.1651. If /tot/-^iayor = 1 nA, then we 
expect Vh ~ 0.16nV, which might be experimentally accessible. Under the same conditions the conventional Bernoulli 
signal is considerably smaller, F ~ 1 pV. It is worth noting that the capacitor technique in measuring the transverse 
potential change does not require a Hall bar geometry, but a single surface is sufficient. 

Unfortunately, besides the high sensitivity necessary in this kind of measurements a further problem has to be taken 
into account. This is the formation of domains of degenerate superconducting states with opposite chiralities. Such 
domain formation is very likely to occur when a system enters the superconducting state. The two domains with 
opposite chiralities yield opposite sign of the transverse voltage so that the net effect might be diminished. Since the 
spontaneous Hall voltage is a surface effect, the number of domain walls intersecting the surface matters. It would be 
necessary to establish an experimental technique to realize a single domain phase, for example, by cooling in a weak 
magnetic field. 
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VII. NON-CHIRAL TIME-REVERSAL SYMMETRY BREAKING STATES 



We now consider the possibility of a spontaneous Hall effect in non-chiral time reversal symmetry breaking super- 
conductors. While for these superconducting states there is no anomalous CS-like coupling between scalar and vector 
potential, there are still spontaneous surface currents for certain orientations of the sample boundaries, despite the 
fact that the Cooper pairs do not have a net angular momentum. These surface currents can be associated with 
Andreev bound states. Thus we would expect at least to find a contribution to the SHE due to the Bernoulli force. 
We consider one of the well-known examples of a non-chiral superconducting state in a quasi-two-dimensional system 
that violates time-reversal symmetry, the d + is-wave state. We assume that the d-wave state has-the d^jiy-symmetry, 
for which a spontaneous surface current runs along surfaces normal to the [100] or [OlOj-direction.L ^ It is possible here 
to express the spontaneous current again as a result of the order parameter texture analogous to (|6.3[) : 

vo = f = ^ f j^'^ ! "^i^/H ) > (7-1) 

where rjs and rjd denote the order parameters of the s-wave and d-wave component respectively. The currents can be 
of similar magnitude as for the chiral p-wave state and consequently the size of the SHE is comparable. There are, 
of course, some differences from the former case. Since the two pairing components are not degenerate in general, 
various additional parameters may appear in the discussion. We have performed a BdG calculation for a specific set 
of parameters to verify the expectation of the above argument. 

Figure ^ shows the data obtained from the self-consistent solution of the BdG equation for the d^y + zs-wave state 
at zero temperature. Here we present the results for the case where a finite net current is running in the system. 
Again the order parameter shows strong variation at the surface, whereby the d-wave component is suppressed and the 
s-wave component enhanced [Fig. §(a)]. Looking at the quasiparticle spectrum in Fig. §(b), we see obvious differences 
between the chiral and the non-chiral cases. In both cases there are chiral Andreev bound states below the ordinary 
continuous spectra of scattering states. In the non-chiral case, however, there is no gapless edge mode, in accordance 
with the expectation from the index theorem. The electromagnetic properties in the Hall bar geometry, shown in 
Figs. |(c) and (d), are very similar to the chiral case. Also, the Hall resistance Rh for the parameters indicated in the 
figure caption is of similar magnitude. We thus conclude that for certain surface orientations the measurement of the 
transverse voltage in a non-chiral state gives a qualitatively identical result to a chiral state. Hence, experiments of 
this kind on the [110] surface of high-temperature superconductors, where a low-temperature time-reversal symmetry 
breaking phase may be present, would not be able to give decisive results as to which state is realized, the non-chiral 
dx2_y2 + is-wave state or the chiral dx^-w^ + ida:„-wave state. 



VIII. CONCLUSIONS 



We have analyzed in detail the spontaneous Hall effect in time reversal symmetry breaking quasi-two-dimensional 
superconductor of chiral and non-chiral nature. There are two contributions to the SHE. One is connected with 
the Bernoulli (or Lorentz) force due to the presence of spontaneous surface currents. The other originates from the 
presence of an orbital angular momentum of Cooper pairs. While the former contribution appears in both types of 
superconducting states, there is no angular momentum in the non-chiral case. We have shown in our phenomenological 
treatment that the angular momentum gives rise to a Chern-Simons-like term in the Lagrangiau. determining the 
electromagnetism of the superconductor, similar to derivations based on topological arguments. Q'E3 Although it was 
suggested that the Hall response would be, at least, close to a universal value, our analysis shows that the actual 
measurement of the Hall voltage gives a considerably smaller non-universal value. Nevertheless, the comparison of 
the SHE with the ordinary Hall effect reveals the presence of an intrinsic effective magnetic field corresponding to the 
density of Cooper pair magnetic moments. 

The effect depends on the presence of spontaneous surface currents. In our analysis we have restricted ourselves to 
the case of a perfect, specularly scattering surface. Rough surfaces with diffuse scattering would reduce the spontaneous 
currents and the Hall voltage. Moreover, domain formation constitutes another obstacle to the measurement of the 
SHE, because different domains give contributions of opposite signs. Nevertheless, we believe that experimental 
techniques available at present are sufficiently accurate to observe the spontaneous Hall effect. 
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FIG. 1. Self-consistent solution of the Bogoliubov-de Gennes equation of the px + ipj,-wave state at T = 0. The set of 
parameters are chosen as = hvp/Ao, kp^o ~ 16, uic = 8Ao, k — A/^o — 1, £ — 0.25^o and Lx ~ Ly = 20^0. (a) Order 
parameter scaled by Ao, which is the magnitude of the order parameter in the bulk region, (b) Energy eigenvalues obtained 
for sine basis functions. L and R denote the surface bound states localized near x = and x = Lx, respectively. States with 
\En\ > Ao are extended. Unit of ky is 2-K/Ly. (c) Dimensionless vector potential ay = {e^o/fic)Ay, magnetic field bz = ^odxtty 
and current density jy = —^odxbz- (d) Dimensionless scalar potential ao = eAo/Ao, electric field Cx = —^odxao and charge 
density p = ^odxCx- 
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FIG. 2. Self-consistent solution of the Bogoliubov-de Gennes equation at T = when current J = SeAo/Tt/cF^o is externally 
supplied. Set of parameters are the same as used in Fig. 0. (a) Dimensionless vector potential, magnetic field and current 
density, (b) Dimensionless scalar potential, electric field and charge density. 
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FIG. 3. (a) K dependence of the Hall resistance. Within the same kp^o, the three lines from the top correspond to ^ = O.l^o, 
£ = 0.25^0 and £ = 0.5^0, respectively, (b) £ dependence of the Hall resistance. Within the same /cfCo, three lines from the top 
correspond to k — 1, k — 1.5 and k = 2, respectively. In both figures the Hall resistance R is scaled by -Ro = h/2e^. 
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FIG. 4. Self-consistent solution of the Bogoliubov-de Gennes equation for the [100] surface in the d^y + is-wa.ve state at 
T = 0. The set of parameters are the same as used in Fig. ^. (a) Order parameter scaled by Ao, which is the magnitude of the 
order parameter for the d^^y-wave in the bulk region. Ad and As are the magnitude of the order parameter of the dxy-wave and 
s-wave, respectively, (b) Energy eigenvalues obtained for sine basis functions. L and R denote the surface bound states localized 
near left and right side surfaces, respectively, (c) Dimensionless vector potential Uy = [eS^o/hcjAy, magnetic field = S,odxay 
and current density jy = —S^odxbz- Total current is / = beAo/hkpS^o- (d) Dimensionless scalar potential ao — eAo/Ao, electric 
field — —(,odxao and charge density p — ^odx^x- 
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